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Abstract 

We study a model of radiating gases that describes the interaction of an inviscid gas with 
photons. We show the existence of smooth traveling waves called 'shock profiles', when the 
strength of the shock is small. Moreover, we prove that the regularity of the traveling wave 
increases when the strength of the shock tends to zero. 

1 Introduction and main results 

We are interested in a system of PDEs describing astrophysical flows, where a gas interacts 
with radiation through energy exchanges. Similar questions arise in the modeling of reentry 
problems, or high temperature combustion phenomena. The gas is described by its density 
p > 0, its bulk velocity n G R, and its specific total energy E = e + where e stands for 

the specific internal energy. (Our analysis is restricted to a one-dimensional framework, but this 
is not a loss of generality, as shown below.) We consider a situation where the gas is not in 
thermodynamical equilibrium with the radiations, which are thus described by their own energy 
n. The evolution of the gas flow is governed by the system: 

dtp + dx{pu) = 0, 

dt{pu) + dx{pu'' + P) = Q, (1) 
dt{pE) + dx{pEu + Pu) = n-9^, 

where the right-hand side in the last equation accounts for energy exchanges with the radiations, 
P being the pressure of the gas, and 9 its temperature. Throughout the paper, we always assume 
that the gas obeys the perfect gas pressure law: 

p = Rp9 = {-f -l)pe, (2) 

where R is the perfect gas constant, and 7 > 1 is the ratio of the specific heats at constant 
pressure, and volume. This assumption yields many algebraic simplifications, but we believe 
that our results still hold for a general pressure law satisfying the usual requirements of thermo- 
dynamics. System is completed by considering that radiations are described by a stationary 
diffusion regime that reads: 

-dxxn = 9'^ -n. (3) 

We detail in Appendix EI how the system (Q), Q can be formally derived by asymptotics 
arguments, starting from a more complete system involving a kinetic equation for the specific 
intensity of radiation. 



As a matter of fact, the operator (1 — dxx) can be explicitly inverted, and Q can be recast 
as a convolution: 

n{t,x) = \ [ e-^^-y^e{t,y)'dy. (4) 



(5) 



Let us introduce the quantity: 

q{t, x) := -dxU {t, x) = l e-l'^-^l sgn(x - y) e{t, yf dy , 
where sgn is the sign function: 



sgn^xj 



1 if X > 0, 
if X = 0, 
-1 ifx<0. 



The quantity q can be interpreted as the radiative heat flux. Then, we can rewrite Q as 
follows: 

' dtp + dx{pu) = 0, 



dt{pu)+dx{pu' + P) = {), 
dt{pE)+dx{pEu + Pu + q) =0, 



(6) 



with q given by ©• Recall that E = e + /2, and P is given by ((3). 

In this paper, we address the question of the influence of the energy exchanges on the 
structure of shock waves. More precisely, let us consider given states at infinity (/o±, Ui, e±), 
and let us asume that: 



{p,u,e){t,x) 



(yO_,M_,e_) ifx<(Tt, 
M+, e+) if X > at, 



(7) 



is a shock wave, with speed a, solution to the standard Euler equations (that is, system © with 
q = Q). We refer to |Lax73l ISer99l ISmo94j for a detailed study of shock waves for the Euler 
equations. The question we ask is the following: does there exist a traveling wave [p, u, e){x — at) 
solution to (jnj, with q given by ©, that satisfies the asymptotic conditions: 



lim (/9,n,e)(^) 

5— »±oo 



{p±,u±,e±) 



(8) 



In other words, we are concerned with the existence of a shock profile, and a natural expectation 
(at least for shocks of small amplitude) is that the step shock ((7j) is smoothed into a continuous 
profile, due to the dissipation introduced by ©• The analogous problem for the compressible 
Navier-Stokes system has been treated a long time ago, see |Gil51j . without any smallness 
assumption on the shock wave. Concerning radiative transfer, a formal analysis of shock profiles 
has been performed in |Hea63j . together with rough numerical simulations. (We refer also to 
|Zel66[ IMih84j for the physical background.) The main purpose of this work is to make the 
analysis of |Hea63j rigorous. Since we are only concerned in this paper with the existence of 
shock profiles, and not with their stability, the problem is purely one-dimensional (due to the 
Galilean invariance of the Euler equations). This is why we have directly restricted to the 
one-dimensional case. However, the formal derivation of Appendix ^ is made in several space 
dimensions. 

Before stating our main results, let us mention that a simplified version of © has been 
introduced, and studied in |ST92j and later in |Kaw99bl IKaw99a] . This 'baby-model' consists 
in a Burgers type equation: 

dtu + dx(^Y) = -^2:9, 

coupled to the diffusion equation: 



-dxxQ + q = 



—dxU . 



These two equations can be seen as a scalar version of ©, © since they can be recast as: 

5tn + 9,(y) =i^n-n, (9) 
where K is the integral operator already arising in 0: 

Ku{t,x) = \ [ e-\''-y^u{t,y)dy. 

The thorough study of © has motivated a lot of works; we mention in particular jNisOOl ILat03l 
ILhiOH ISerOSj . Clearly Q can be seen as a prototype for discussing nevertheless, 
replacing ©, and © by Q has two important consequences: the equation becomes scalar, and 
the 'diffusion' K — \ applies to the unique unknown (while in ©, the 'diffusion' appears through 
the radiative heat flux q only in the third equation). Our work is a first attempt to extend the 
known results for © to the more physical model ©, ®- 

Let us now state our main results. The first result deals with the existence of smooth shock 
profiles when the strength of the shock is small: 

Theorem 1. Let 7 satisfy 

\/7 + 1 

1< 7 < 2.215, 

V7-1 

and let {p-,u-,e-) be fixed. Then there exists a positive constant 6 (that depends on {p-,u-,e-), 
and j) such that, for all state (/9+,-u+,e+) verifying: 

• \\{p+,u+,e+) - {p-,u-,e-^)\\ < 6, 

• the function ^ is a shock wave, with speed a, for the (standard) Euler equations, 

then there exists a traveling wave {p,u,e){x — at) solution to ©, ©, 

As in the study of the 'baby- model' @, the existence of a smooth shock profile is linked to 
a smallness assumption on the shock strength, see |Kaw99b] . Here the smallness parameter 6 
may depend on the state (/9_,ti_,e_), while for (jUj), the smallness parameter is uniform (and 
even explicit!). 

The restriction on the adiabatic constant 7 might be unnecessary, but it simplifies the proof, 
and it covers the main physical cases 1 < 7 < 2. 

Our second result is also in the spirit of |Kaw99b] . and deals with the smoothness of the 
shock profile constructed in the previous Theorem: 

Theorem 2. Let 7 satisfy 

\/7+ 1 

1< 7 < ~ 2.215, 

V7-1 

and let (p-,u-, e_) be fixed. Then there exists a decreasing sequence of positive numbers {6n)nm 
(the sequence depends on (p_,ii_,e_), and 7J such that, for all n E N, and for all state 
(p+,n+,e+) verifying: 

• \\iP+,u+,e+) - {p^,u^,e^)\\ < dn, 

• the function Q zs a shock wave, with speed a, for the (standard) Euler equations, 

then there exists a C"^^ traveling wave {p,u,e){x — at) solution to @, (©• 

To a large extent, our analysis follows the arguments of |Hea63j . |ST92j and |Kaw99b| . The 
proof of Theorem n is presented in Section [21 while Section is devoted to the proof of Theorem 
121 The investigation of strong shocks, as well as stability issues will be addressed in a forthcoming 
work. 



2 Existence of smooth shock profiles 



In this section, we prove Theorem^ We first recall some basic facts on shock waves for the Euler 
equations. Then, we make some transformations on the traveling wave equation. Eventually, 
we prove Theorem ^ by using an auxiliary system of Ordinary Differential Equations, that is 
introduced and studied in the last paragraph of this section. 



2.1 Shock wave solutions to the Euler equations 

In this paragraph, we recall some basic facts about the (entropic) shock wave solutions to the 
Euler equations: 

' dtp + dx{pu) = 0, 
< dt{pu) + dx{pu^ + P)=^, 
dt{pE) + dx{pEu + Pu) = Q, 

where P, and E are given as in the introduction. We refer to |Lax73[ ISer99l ISmo94j for all the 

details, and omit the calculations. In all what follows, we only consider shock waves that satisfy 
Lax shock inequalities. We shall thus speak of 1-shock waves, or 3-shock waves. 

We consider a fixed 'left' state (/>_,«_, e_). Then the 'right' states (/9+,u+,e+) such that 
{p±-,u±,e±) define a 1-shock wave, with speed a, is a half-curve initiating at (p_,u_,e_). In- 
troducing the notation v± = u± — a, the Rankine-Hugoniot jump conditions can be rewritten 
as: 

v+= p^v^ =: j , 

p- v'i =: J L.1 



p+ vl + P+= p- vl + P =: j Ci 



(10) 



P+ 



/9„ V- 



+ —)+P_v^ 



--■■jQ 



2 ■ 



Observe that V- does not only depend on the 'left' state (/9_,M_,e_), but also on (/)+, e+), 
because V- is defined with the help of the shock speed a. Consequently, the constants j, Ci, 
and C2 depend on both (p_,u_,e_), and (p+,'u+,e+). 
For 1-shocks, that is when the inequalities: 



< a < Ua 



a < 



are satisfied (c denotes the sound speed), all quantities j, Ci, and C2 are positive. Moreover, 
when the strength of the shock tends to zero, that is when (p-|_,n+,e+) tends to (/0_,ti_,e_), 
one has 

/ u^-c- \ 
p- c_ 
+ (7 - l)e_/c- 



3 

Ci 

\C2j 



\ 



(11) 



Consequently, when the strength of the shock is small, all quantities j, Ci, C2 are bounded away 
from zero. Recall also that 1-shocks are compressive, in the sense that p^ > p-. This inequality 
immediately implies that < f + < V-. Eventually, the strength of the shock tends to zero if, 
and only if tends to n_ (in that case, we also have p+ — > p_, and 64. e_ because of the 
Rankine-Hugoniot jump conditions). 

In all what follows, we limit our discussion to the case of 1-shocks for simplicity, but the 
extension to 3-shocks is immediate. 



2.2 Reduction of the traveling wave equation 

In this paragraph, we derive, and transform the equation satisfied by traveling wave solutions 
to (0). A traveling wave solution to ©, © with speed a is a solution (p, u, e){x — at). For 
such solutions, the radiative heat flux q also depends on the sole variable x — at: 

q{x -at) = l [ e-\^-^'-y\ sgn(x - at - y) 9{yf dy , 



and © reads: 

'{p{u-a)y = 0, 
< (/9 n (n - cr) + (7 - 1) /9 e)' = , 
Xpie + ^){u-a) + {j-l)peu + qy = 0, 

where ' denotes differentiation with respect to the variable ^ = x — at. Introducing the new 
unknown v = u — a, the above system is easily seen to be equivalent to: 



iipvy = o, 

ipv^ + i^-l)pey = 0, 
[ipvie + ^) + {j-l)pve + qy 



(12) 



0. 



Since we are looking for a shock profile, the traveling wave solution should also satisfy: 

lim {p,v,e) = {p±,v±,e±) , (13) 

where v± = u± — a, and {p±,u±, e±) defines a 1-shock wave with speed a for the Euler equations. 
Notice that the quantity a = — n_|/2, that measures the strength of the shock, is invariant 
with respect to our change of velocity, that is a = — U-\/2 = \v+ — v-\/2. Recall also that 
for 1-shocks, there holds u _ > 1'+ > 0. 
Observing that we have: 



1 r+00 

9(6 = 77 / e-y{e{C-y)'-9{^ + y)^)dy, 
^ Jo 



we conclude that q tends to zero at ±00 by Lebesgue's Theorem (because 9 is necessarily bounded 
since it has finite limits at ±00). We can thus integrate the system H12|) - H13|) once, and ((T^ 
reads equivalently: 

'ipvm=j, 

< {pv^ + {j-l)pe){0=jCi, (14) 
^{pv{e + ^) + {j-l)pve + qm=jC2, 

where the constants j, Ci and C2 are given by the Rankine-Hugoniot conditions (|l()j) . For small 
shocks, the positive constants j, Ci, C2 have the asymptotic behavior ((TT|) . 
Prom the two first equations of (|T1|) . we derive the relations: 



J 



The third equation of 1)141) thus reduces to: 



^^^^2(7 



{Ci-v{0)v{0 
7-1 



1)C2 2(7-1) 



7 + 1 ■ 7 + 1 J(7 + 1) 

Using the equation of state Q, as well as the second equation of (fTl|) . we get: 

= R = R • 

Consequently, (P3|) can be recast as an integral equation with a single unknown function v. 

(7-1) 



(15) 



7 + 1 



7 + 1 



j (7 + 1) i?4 



sg-a{i-y)v{y)\Ci-v{y)fdy. (16) 



We are searching for a solution v to (|16() . that satisfies the asymptotic conditions v{^) v±, as 
i ±00. 



Remark 1. If we find a solution v to (|16() that does not vanish, and that satisfies v{£,) — > v± 
as ^ ^ ±00, then we obtain a shock profile {p, u, e) by simply setting: 

In particular, if v{^) G for all ^, then v does not vanish. 

Remark 2. Since the heat flux q vanishes at ibcxD, ()16() can be also rewritten as: 

{viO - v-){v{0 - v^) = / e-\^-y\ sgn(^ - y) v{yf {C, - v{y))^ dy . 

J (7 + 1) W' 

We are going to rewrite 1)16^ as a second order differential equation, that will be easier to 
study than the integral equation (|16|) . Indeed, assuming that is a C'^ function of ^, and differ- 
entiating twice (fTC)) with respect to we get (see |Hea63j for the details of the computations): 

- ^" + (^')' - ^ (^1 - - 2v) v'-\{v- V.) {v-v+) = 0. (17) 

Conversely, if is a solution to (|17|) that satisfies v{^) — > v-t as ^ — > ±00, then v is also a 
solution to H16() . If in addition v takes its values in the interval then we can construct 

a shock profile, and thus prove Theorem ^ 

The differential equation (|17|) can be simplified by introducing the new unknown function 
V = V — {v- + v+)/2, and by rewriting the second order differential equation as a first order 
system: 



v' = w , 

vw = —w—f[v)w-\ 



(18) 



where / is the following polynomial function: 

' jR'{j + l)\l + l J \ 7 + V V 1 + 1 J ^ ' 

We recall that a = — v^\/2, and that a measures the strength of the shock. 

Remark 3. The asymptotic behavior (|1H) ofj, Ci, and C2 shows that when the strength of the 
shock tends to zero (a — > 0"''^, the limit of f{0) is given by: 

, 473(7-1)2 (c_ + (7-1) ^)^ 
RH-f + lf p-c- 

Since V-^. < V- for a 1-shock, we are searching for a solution to (fTH)) that is defined on all M, 
and that satisfies: 

lim ({},-(/;) (^) = (a, 0) , lim (u, u;)(^) = (-a, 0) . (20) 

To prove Theorem we are thus reduced to showing the existence of a heteroclinic orbit 
for (|18|) that connects the stationary solutions (±0,0). Due to the previous transformation 
V = V — {v- +t'+)/2, if V takes its values in [—a, a], then v = v + (f_ +v+)/2 will take its values 
in the interval and therefore will not vanish. 

Remark 4. The system (|18|) is 'singular' at v = 0. Nevertheless, we are searching for a smooth 
solution connecting (±0,0), so that v vanishes in at least one point. Because w' = v" should 
also have a limit at this point, a C"^ shock profile can exist only if the equation: 

u;2 + /(0)ti;±y = 0, 

has real roots. The corresponding discriminant condition turns out to be much less simple than 
the one found in lKaw99b}l for the 'baby model' ©. (In particular, /(O) depends on the shock 
through the constants j, and Ci). This is a first 'nonexplicit' restriction on the shock strength 
to derive the existence of a smooth shock profile. 



Due to the singular nature of the system ()18() at -u = 0, it is more convenient to work on an 
auxihary system of ODEs, where the singularity has been eliminated (at least formally) thanks 
to a change of variables. This procedure was already used in |Kaw99bj . In the next paragraph, 
we shall introduce this auxiliary system, and complete the proof of Theorem^ 



2.3 Existence of a heteroclinic orbit 

We begin with a result on an auxiliary system of ODEs, where the singularity at v = has been 
eliminated: 

Proposition 1. Assume that 7 satisfies 1 < 7 < (\/7 + l)/(\/7 — 1), md consider the following 
system of ODEs: 

V' = VW , 

There exists a positive constant oq, that depends only on (/>_,u_,e_), and 7 such that if the 
shock strength a satisfies a G (0,ao], the following properties hold: 



• /(0)2 - 2a2 > 0, and we define wq:={- /(O) + v^/(0)2 - 2a2)/2 < 0. 

• There exists a solution (V^, to (|21j) that is defined on all M, and that satisfies 

hm {V,,W^){ii) = (a,0), hm {V,,WM = i^.w^) ■ 

»?— 00 77— >+oo 

Furthermore, VJ, is decreasing, and the convergence of to as 7/ — > +00 is exponential. 

• There exists a solution (V{(, Wjj) to (|21j) that is defined on all M, and that satisfies 

hm (yp,VFj)(7?) = (-a,0), hm (T^, TF„)(r?) = (0, u;o) • 

Furthermore, Vj is increasing, and the convergence ofV^ to as rj ^ +00 is exponential. 

Assuming that the result of Proposition^holds, the existence of a heteroclinic orbit for H18|l 
connecting (=ba, 0) can be derived by following the analysis of |ST92| lKaw99bj . We briefly recall 
the method. Using the solution iy^, W^^), we introduce the change of variable: 

r+co 

Since V\, tends to exponentially as r] tends to +00, H[, is well-defined, and it is an increasing 
C°° diffeomorphism from M to (—00, 0). Then (w^; •= (Hj ^b) ° '='b~^ ^ solution to ((TH|l 
on the interval (— oo,0), that satisfies: 

lim (4,ri;J(^) = (a,0) , lim (%, ii;[,)(^) = (0, u;o) • 

Similarly, with the help of the solution (Vf|,Tyn) we can construct a C°°, decreasing diffeomor- 
phism H|j from M to (0, +cxd), and a C°° solution {v^,w^) to (|T8|) on the interval (0,+oo). This 
solution {v^,w^) connects {0,wo) and (— a, 0), as ^ varies from 0+ to +00. Let us now 'glue' the 
solutions {v\f,w\,), and {v^,w^), by defining: 

{v,w){C) := < (22) 

and extend the functions v, and at by setting (v,w){0) = (0, wq). In this way, v, and w are 
continuous on M, and C°° on M \ {0}. It remains to show that v G C^(M), that (u, w) solves ()18() 
on M, and that v takes its values in {—a, a). 




Observe first of all that u is a decreasing function, because of the monotonicity properties of 
, Vjj , H[, , Sjj . Using the asymptotic behavior of V[,, Vj at — oo, we get that v{S,) G (—a, a) for all 
C G M. 

Let us now note that the above construction of {v,w) shows that {vjw) is a solution to H18|) 
on M \ {0}. In particular, v'{S,) = w{(,) if 7^ 0. Moreover, w is continuous on R, so {) G C^(M), 
and v'{0) = w{0) = wq. To prove that v £ C^(M), it is sufficient to show that w G C^(M), which 
is equivalent to showing that w' has a limit at (because we already know that w is C°° on 
M \ {0}). To prove that w' has a limit at 0, we are going to study the asymptotic behavior of 
(14, VFi,), and (VJj, Wj) at +00. More precisely, let us denote U{V,W) the vector field associated 
to the ODE (EH): 

/ VW 

where / is given by H19|) . The Jacobian matrix of U at (0,wo) is: 

Wo 
-f'{0)wo -2wo-f{0) 

For a sufficiently small, one checks that X^2^ < A^*''* < (see Proposition ^ for the definition of 
Wo)- The eigenvectors corresponding to the eigenvalues A^'^^ and X^^ are: 

The standard theory of autonomous ODEs, see e.g. |Pon75j . shows that there are exactly two 
solutions to (PT|) that tend to (0,u;o) as rj tends to +00, and that are tangent to the straight 
line {0,wo) + M.e'^K Moreover, all the other solutions to ^IT\i that tend to {0,wo) as rj tends to 
+00 are tangent to the straight line {0,wo) + We^\ Now, it is rather simple to see that the two 
solutions to (|2H) that tend to {0,wo) as r] tends to +00, and that are tangent to the straight 
line (0, Wq) + 1^62'''', satisfy V = 0, and: 

2 

Because the solutions (H)l^l))> ^-iid (V{|,Ty|j) given by Proposition ^ cannot satisfy V\, = 0, and 
V(| = 0, we can conclude that the solutions (V[,, W;,), and (Vj, Wjj) are tangent to {0,wo) +Me^''^ 
as r] tends to +00. In particular, this yields: 

EM = li^ = -/^(Q)-o (24) 

v^+^V^{v) v^+ooV^iT]) f{o)+3wo' ^ ' 

(A quick verification shows that /(O) + 3 > for small enough a.) Prom the construction of 
the solutions (wi,, it^i,), and (wjj, ifjj), we get: 

lim wliO = lim w[{0 - "^'^^^'^o 



W' = -W"^ - f{fd)W -'^ . 



€-0- i^o+ /(0) + 3w;o 

As a consequence, when a is small enough, w G C^(R), and therefore "0 G C^(M). Moreover, 
{v,w) solves ()18() on M \ {0}, so by continuity, it solves (|18() on R. This completes the proof of 
Theorem n provided that the result of Proposition^ holds. 

2.4 Proof of Proposition [T] 

In this paragraph, we prove Proposition^ which will complete the proof of Theorem^ At first, 
we define the set: 

P={{V, W)\V G [-a, a],W^ + f{V) W - ^^^^^ = o} , 




P9 



Figure 1: The set P = P1LIP2. 

so that the points (±a, 0) belong to P. The following Lemma gives a description of P for a > 
small enough. We refer to figure ^ for a schematic picture. 

Lemma 1. Assume that 1 < 7 < (\/7 + l)/(\/7 — !)• T/ien t/iere exists a constant ao > 0, that 
only depends on {p-,U-, e_) and 7 suc/i that if the shock strength a satisfies a G (0, ao], we have 
the following results: 

• For all V G [-a, a], f{Vf + 2 (F^ _ g^) > 0. We can thus define 

-fiV) + Vf{Vr + 2{V^-a^) 



\/Ve[-a,a], Wi{V) 

W2{V) 



-f{V)-Vf{Vr + 2{V^-a^) 



• P = PiU P2, where Pi and P2 are two curves defined by 

Pi = {{V,Wi{V))\V G [-a,a]} , P2 = {{V,W2{V))\V G [-a,a]} , 
so that the points (iba, 0), and {0,wo) belong to Pi. 

• There exists a unique point V G (— a,0) such that Wi is increasing on the interval [V,a], 
and Wi is decreasing on the interval [—a, V]. 

• For all V G [F,0], one has W2{V) < Wi(F). 
Proof. Let us first define a function A by setting: 

A{V) := f{Vf + 2{V^ -a^). 
Using for a small enough, we have: 

Ci jCi h-l)Ci 

a>K>0, — - a > K > , — - 2a > K > , 



7 + 1 



7 + 1 



7 + 1 



where At is a positive constant that only depends on (/9_,u_,e_) and 7. Moreover, Ullj) also 
shows that j > k for a S (0, ao], up to restricting k. Consequently, there exist oq > 0, and k > 
such that for a £ (0, oq], we have f{V) > k, and A(y) > k for all V S [—a, a]. This directly 
shows that the set P is the union of the two curves Pi, and P2. It is rather clear from the 
definition of Pi that (iba,0), and {0,wq) belong to Pi (recall that wq is defined in Proposition 
d). Observe also that W2{V) < Wi{V) < for F e [-a, a], and Wi{V) < if F G {-a, a). 
The functions Wi, and W2 are C°° on [—a, a]. Moreover, we compute the relation: 



VFG[-a,a], y/A{V)W[{V) = V-Wi{V)f'{V), (25) 
and from (|T!H) . we also compute 



iP4 (^ + 1) 



Ci 7Ci\V. (7-1)Ci Ci\ / . (7-l)Ci Ci \ 



7 + 1 / V 7 + ly V 7 + 1 \/7y V 7 + 1 V7 

As we have done for /, and A, a careful analysis shows that for 1 < 7 < {^/7 + 1)/(\/7 — 1), and 
for a small enough, one has f'{V) > k > for all y G [—a, a], because each term in the product 
(|26|) is positive. Using this information in (|25|) . we can already conclude that Wi is increasing 
on the interval [0, a] (see figure P). Moreover, the relation (^5]) also shows that W'^(O) > 0, and 
W'i{—a) < 0. Consequently, there exists some V G (— a, 0) such that W'^(F) = 0. Let us prove 
that V is the only zero of W'^. We claim that it is sufficient to show the following property: 

W[{V) = ^ W'liV) > . (27) 

Indeed, if the property H27() holds true, then any point where W'^ vanishes is a local strict 
minimum. If there existed two such local strict minima —a < Vi < V2 < a, then Wi would 
admit a local maximum V3 G (^^1,^2)) which is obviously impossible. Therefore let us prove 
that the property (P7|) holds true. 

Differentiating ()25|) with respect to V, we obtain that if W']^(y) = 0, then 

^A(F)w'i'(F) = 1 - /"(F) Wi(F) . 

Observing that 

|/"(F)Wi(F)| < C|Wi(F)| <c^--^< ^ , 

for suitable positive constants C, and k (that are independent of a G (0,ao]), we can conclude 
that W'/(y) > 0, provided that a is small enough. This completes the proof that W'^^ has a 
unique zero V G (—a, 0), and therefore Wi is decreasing on [—a, V], and is increasing on [F, a]. 
For the last point of the lemma, we use the relation: 

W'i(T/) + W2(y) = -f'{V) < 0. 

Because W{{V) > O_for F G JF,0], we get W^(y) < for y G [F,0]. Thus for V G [F,0], we 
have W2(l^) < W2(T^) < Wi(y), and the proof of the Lemma is complete. □ 

Using Lemma ^ we are going to prove Proposition ^ The analysis follows |Gil51j . 

As we have already seen in the preceeding paragraph, the point {wq,^) is a stable node of 
H21|) . We now study the nature of the equilibrium points (±a, 0). Recall that the vector field 
associated to (|^T|) is denoted U, see ((221) • The Jacobian matrix of U at (a, 0) is: 

a 

a f{a] 




Figure 2: The compact set Ki. 



so it has exactly one negative eigenvalue ni, and one positive eigenvalue iJ,2 (the equilibrium 
point (a, 0) is a saddle point): 



-/(a)-x//(«)^ + 4a^ -/(«) + V/(a)^ + 4a^ 
W = ^ , /^2 = . 

An eigenvector associated to fj,2, and is r2 = (0,^2)- Moreover, using the relation (^5]) . we can 
check that for a small enough, the following inequality holds: 

0<^<-f- = W;(a), (28) 

where the function Wi is defined in Lemma ^ Let us now define a compact set Ki by: 

Ki := |(y,VF) G [0,a] x R| Wi(y) < < o} , 

Then the inequalities H28|) show that for s < small enough, the point (a,0) + sr2 belongs to 
the interior of Ki. We refer to figure I^J for a detailed picture of the situation. 

From the standard theory of autonomous ODEs, see e.g. |Pon75) . we know that there exists 
a maximal solution (F^,Wi,) to 1)21(1 that tends to the saddle point (a,0) as rj tends to —00, 
and that is tangent to the half-straight line (a, 0) + M~ r2- This solution is defined on an open 
interval (— oo,ry*) (with possibly ??* = -|-cxd). For large negative rj, the preceeding analysis shows 
that (y\>,W\,){r]) belongs to the interior of Ki. Moreover, (H)^b) cannot reach the boundary 
of Ki. Indeed V\, cannot identically vanish so (H;^b)(^) 9Ki n {V = 0}. Similarly, we have 
(y\,,W\,){r]) 7^ (a, 0). Eventually, on the set: 



{(y,o)|y G (o,a)}u{(y,Wi(y))|y e (0,a)} 



the vector field U is not zero, and is directed towards the interior of Ki. Therefore the solution 
(V^, W(,) cannot reach dKi, so it takes its values in the compact set Ki. The maximal solution 
(yi,, W[,) is thus defined on M. It cannot reach the boundary of Ki, so W[, takes negative values, 
which means that V[, is decreasing (because = V\,Wi,). Because (Vb,W[,) takes values in the 
interior of Ki, the function W\, is also decreasing. This shows that (y\,,W\,){r]) has a limit as r] 
tends to -|-cxd, and this limit is necessarily be a stationary solution of ()21j) . The only possibility is 
that {V\,, W\,){r]) tends to {0,wo) as r/ tends to +00. The convergence is necessarily exponential, 
because the Jacobian matrix of U at {0,wo) has two negative eigenvalues, see e.g. |Pon75j . 
To construct the other solution (Vj, W^), we argue similarly by defining a compact set K2: 

K2 := \^{V,W) G [-a,V] x M|Wi(y) < < o} U |(F,VF) G [F, 0] x M|Wi(F) < < o} , 

see figure 131 The Jacobian matrix of U at (— a,0) has one negative eigenvalue vi, and one 
positive eigenvalue 1^2, with: 




Figure 3: The compact set K2- 



An eigenvector associated to the eigenvalue 1^2 is R2 = (~Q^)^2)- As was done earlier, we check 
that the inequalities: 

/(-a) -a 

hold true. Therefore, one can reproduce the above analysis, and show that there exists a solution 
(Vjj, VFjj) to (PT|) that takes its values in K2 (and is thus defined on M), and that tends to (—a, 0) 
at —00. Moreover, (Vij, W^) can not reach the boundary of K2, so is increasing. It only remains 
to study the monotonicity of Wf^. This is slightly more complicated than for W\,. Observe that 
K2 is the union of the sets: 

Kl := G [-a,0] x M| Wi(y) < < o} , 

■= ^{V,W) G [F,0] X M| Wi(F) <w < Wi(y)} . 

When (Vj, W(j) takes its values in the interior of K2, the function VFj is decreasing (this is the 
case for large negative r]). At the opposite, when {Vij, W^) takes its values in the interior of K2, 
the function VFj is increasing, because thanks to Lemma we have: 

= (Wi(ys(r/)) - W^{r^)) (W^iT^) - W2{mr]))) 

> (Wi(ys(r/)) - W^ir])) (Wi(F) - W2(ys(r/))) > 0. 

Moreover, if (VJj, W^){r]o) belongs to the interior of -fC| for some r]o £ M, then (VJ, W^){r]) belongs 
to the interior of -K"! for all rj > r]Q (because it cannot reach the boundary of K2 for r] > ijq). 
Summing up, either (Vjj, W(()(r/) belongs to K2 for all t], and is monotonic on M, either 
(Vj,Wn)(r/) belongs to K2 for all rj greater than some t/q, and is monotonic on [rjQ, +00). 
In any case, the function is monotonic on a neighborhood of +00, and thus has a limit at 
+00. This shows that {Vi^,W^){r]) tends to (0,^0) as t] tends to +00, and the convergence is 
exponential. As a matter of fact, we have seen in the preceeding paragraph that (Vf[,Wn) is 

tangent to the straight line (0, wq) + M e^'^^ as tj tends to +00, so one can check that (Vj, W^){r]) 
belongs to the interior of K2 for large positive r]. This means that VFj is decreasing on some 
interval (—00,7/0), and increasing on [r]o,+oo). The proof of Proposition ^ is now complete. 



3 Additional regularity of shock profiles 

As should be clear from the preceeding section, the key point in the construction of a shock 
profile is Proposition ^ that gives the existence of two heteroclinic orbits for the system ((71]) . To 



prove Theoreml^ we are going to study the behavior of the derivatives of {V\,, W\,), and (VJj, Wj) 
near +00. The proof of Theorem [21 follows from an induction argument. To make the arguments 
clear, we deal with the first case separately. In all what follows, (V^, Wi,), and (Vj,VF|j) are the 
solutions to (PT|) that are defined in Proposition^ and {v,w) denotes the solution to (fTH]) that 
is defined by (|22|) . We have the following: 

Proposition 2. Under the assumptions of PropositionU[ there exists a positive constant ai < 
aQ (that depends on (/9_,u_,e_), and such that for all a € (0,ai], one has w G C^(M), 
V G C3(E), and: 

= WQ + wi v{C) + W2 v{C)^ + o{v{C)^) , as ^ ^ 0, 

for some suitable constants wi,W2 (wq has already been defined in Proposition^. 

Proof. Recall that Vb, and VJt do not vanish on M, so we can introduce some C°° functions W\, 1, 
and Wjj^i that are defined by: 

Wi, = wo + H ^-^,1 , Wi = wo + Vi Wii . 

Substituting in (PT|) shows that (V[,;^b,i)) (^S)^tt>i) solutions to the system: 

'V = V{wo + VWi), 

f(V) - f(0) ^ V (29) 

W{ = -wo ^^ ^ -2VW^-3woW^- f{V) W, + -. 

Moreover, we already know from Proposition^ and (|24() that: 

-no) wo 



hrn (H,t^b,i)(^)= lim (Fj, W^«,i)(r/) = (o. 



We denote Ui{V,Wi) the vector field associated with (|^: 

V{wo + VWi) 

Recall that / is a polynomial function of degree 7, see (fT^ . thus F{V) := {f{V) — f{0))/V is 
a polynomial function of degree 6, and we have F{0) = /'(O), F'{0) = /"(0)/2. Obviously the 
system of ODEs ()29() admits the equilibrium point {0,wi), where: 

-F{0)wo -f'iO)wo 



f{0) + 3wo f{0) + 3wo 

We are now going to study the nature of the equilibrium point (0, wi), and show that for a small 
enough, this equilibrium point is a stable node for ()29|). Then we shall show that w G C^(M), 
and V G C^(M). In the end, we shall derive the asymptotic expansion near ^ = 0. 
Step 1: the Jacobian matrix of Ui at {0,wi) is: 

Wo 

^^-2wf-f{0)w,-^wo 

Using Remak|31 we can conclude that for sufficiently small a, that is a G (0, ai] for some positive 
number ai less than oq, one has A2^'' < a["^^ < 0, that is, /(O) + Awo > 0. Moreover, the 
eigenvectors corresponding to the eigenvalues X^i \ and A2^'' are: 





Consequently, {0,wi) is a stable node of ((^. and there are exactly two solutions to (PUj) that 
tend to {0,wi) as rj tends to +00, and that are tangent to the straight line {0,wi) +Me2^\ All 
the other solutions to (|29|) that tend to (0, wi) as tj tends to +00 are tangent to the straight line 
{0,wi) + Me[^\ As in the preceeding section, we can thus conclude that: 

Step 2: if we let gi denote the second coordinate of the vector field Ui, we have = 
91(H) W^b.i)) and Wjj^ = g'i(Vj, Wjj^i). Differentiating once with respect to rj, and using (|3fl]) . we 
end up with: 

Wl\{r]) W'\i7]) 
lim : = lim ;' = h , (31) 

where the real number ii can be explicitely computed (but its exact expression is of no use). 
Following the analysis of the preceeding section, we define some functions w\, i := W\, i o E^'^, 
and wj^i := Wj^i o H^"^. First of all, (|3()|) yields: 

=^lim 4.(0 = jl|^, (32) 
Observe now that we have the relations: 

and combining with (|31j) . we get: 

^hm (z), ^ J'(0 = ^lirn^ --^ = ^lim^ ' = ^0 , 

lim (-yH 1) (4) = hm = hm — =Iiwq. 

^^0+ " 77-^+03 Vj(r/) »y^+oo T/^'(77) 

Differentiating twice the relations W], = wq + % f^b,!) and tfjj = t(;o + f^tl.i' obtain: 



(33) 



W^^' = v'^ n,l + ^^b ^b,l + (% ■«^b,l)' = "^b ^«b,l + ^«b ^«b 1 + ("f^b n,l)' ' 



w'l = w\ wii + + (% -wj^i)' . 

Using JHSl), and (jSHI), we get w'^{<d') = wj'(0+). Using the definition this shows that 

w G C^(M), and using v' = w, we obtain v € C^(R). 

Step 3: note that we have the following expansions near ^ = 0: 

^(0 = w{0) + w'{0) e + ^ e + oif) , 

^(O = ^(o)c + ^e^ + o(a, 

with w{0) = wq < 0. We can thus combine these expansions, and derive: 

for some appropriate real numbers a, and /?, that we are going to determine. From the relation 
'"^b(0 = ""^0 + 11)1,^1 {^), and using that ifb,i(?) tends to wi as ^ tends to 0~, we first get 
a = wi. Then from H32|) . and from the relation ^^(0~) = wq, we can obtain: 

^^b,i(0 = ^^1 + /I *^(^) + o(^b(0) > as C ^ 0~. 

We thus obtain /3 = 6i/(/(0) + 4i(;o), which yields: 

w^(6 = + w^i 1^(6 + W2 v{0^ + o(t;(C)^) , 

where W2 ■= 6i/(/(0) + 4:Wq). This latter expansion will be generalized to any order in what 
follows. □ 



We now turn to the proof of Theorem |2J More precisely, we are going to prove the following 
result, that is a refined version of Theorem [2J 

Theorem 3. Let the assumptions of Proposition^be satisfied. Then there exists a nonincreasing 
sequence of positive numbers (an)nGN such that, for all integer n, if a G (0,a„], then w G 
C"~''^(M), and v G C"^^(R). Moreover, w admits the following asymptotic expansion near 

where the real numbers wq, . . . , Wn+i are defined by: 



Wo 
Wk 



-f{0) + Vf{0r-2a^ 

Ofc-1 



f{0) + {k + 2)wo 
and the real numbers bo, . . . ,bn are given by: 



for fc = 1, . . . , n + 1, 



bo 



-no) wo, 



1 



b^ = l-2wl-f{0)wi 



no) 



■Wo ■ 



^ — ^ — Wk+i-i - + 1) Wi Wk+i-i , for A; = 2, . . . , n. 



i=l 



i=l 



Proof. The case re = has been proved in the preceeding section, while the case n = 1 is proved 
in Proposition 121 (The reader can check that the definition of wo, wi, W2, bo, and bi coincide 
with our previous notations.) We prove the general case by using an induction with respect to 
re, and we thus assume that the result of Theorem holds up to the order re > 1. We are going 
to construct a„+i so that the conclusion of Theorem |21 holds for a G (0, fln+i]- In particular, the 
real numbers wo, . . . , Wn+i, and bo, . . . ,bn are given as in Theorem |21 and we can already define 
the real number 6^+1 by the formula: 



bn 



+1 



n+2 n+1 
}_^——Wn+2-^ 



Wi Wn+2-i 



i=l 



1=1 



(Observe indeed that this definition only involves wq, . . . ,Wn+i, and not Wn+2-) 

Step 1: because VJ,, and do not vanish, we can introduce some functions and 
W^tt:"+i relations: 



= ^i;0 + W^l H + • • • + + ^b,n+l ' Wi = Wo + WiVi + --- + WnVf + W^^n+l 

Thanks to Taylor's formula, we can write the polynomial function / as: 



n+l 



re! 



f{V)=f{0) + f{0)V + i^V^ + 
where Fn+i is a polynomial function such that: 

Substituting the expression of W\,, and in (PT|) shows (after a tedious computation!) that 
(14) W\,^n+i)j and (Vft, Wn^„+i) are solutions to the following system of ODEs: 



V = V{wo+WiV + ---+WnV'' + Wn+l F"+^) 
VF;+l=9n+l(^,H^n+l), 



(35) 



where the function Qn+i is given by: 

/ n \ n 

gn+l{V, Wn+l) := - (n + 2) Wn+l [Y^WkV^ + 1/"+^ Wn+l - Wn+1 Y,{k + 1) Wk 



\k=0 / k=0 

n 



- Wn+l f{V) - Fn+l{V) y2wkV'' + bn+ \ ^\ +VQn+l{V) , (36) 

and Qn+i is a polynomial function that satisfies: 

(5„+i(0) = bn+i + (n + 4) u;„4.i + /'(O) u;„,+i + - — — — wi + - — — — wq . 

(n + Ij! (n + zj! 

When n = 1, one has (^2 = Q2(0) = (see the above definition for &2)- Using the expansion 
(|34|) ■ which is part of the induction assumption, we also know that: 



lim (H,t^b,n+i)(^) = lim m,Wi^n+i){v) = (0,w;„+i) = (o, 



f{0) + {n + 2)wo 



With the above definitions for gn+i, and Qn+i, we can check that {0,Wn+i) is a stationary 
solution to (|35j) . (Recall that u;„+i is defined as in Theorem 01 by the induction assumption.) 
We can also evaluate the Jacobian matrix of the vector field associated with the system of ODEs 
(|35ll: 

^0 A_AS"+^) ' 



S+i -m -in + 3) wo J \ hn+i a("+') I ■ 

There exists a positive number a„+i < a„, such that for all a G (0,a„+i], one has a("+') < 
^(n+2) ^ equivalents /(O) + (n + 4) -wq > 0. In that case, the eigenvectors corresponding 
to the eigenvalues A^"^"*^^ and A2"^"'^^ are: 

Using the same argument as in the proof of Proposition ^ we can conclude that the solutions 

{V\f,W\,^n+i)j and (H'^tt:"+i) are tangent to the straight line {0,Wn+i) +Me^"^^^ as rj 

tends to +oo. In particular, this yields: 

.i?^oo - ,^^foo -V^^ - f{0) + {n + 4)wo ■ ^^^^ 



Step 2: let us define the function tx^n+i by the formula: 



fvF,,„+ioH-i(e) ife<o, 

Wn+l if C = 0, 

„+ioH-i(e) ife>o. 

With this definition, ivn+i is continuous, and we have the relation: 

W = Wo+WiV-\ h U^n + Wn+l ^""^^ • (38) 

Moreover, using 1)3 7p . we obtain: 

^™ -//c^N = = Wn+2 , (39) 

which yields 'w;^^i(0"'") = w'^_^_i{0~). Therefore, we have Wn+i S C"^(M). Moreover, using 
we can compute: 

^^^n+l V = 9n+l(.V, Wn+l ) , (40) 



so we get w'n^i v G C^(M). 

Step 3: we use an induction argument to show that w S C"'"^^(M) (which will imply im- 
mediately V G C"'^^(IR)). More precisely, we assume that for some k G {0, we have: 

Wn+i G C^+^m , w e C'=+i(M) , v''+^ G C''+\R) . (41) 

We are going to show that this property implies the same property with k replaced by A; + 1. 
(Observe that step 2 above shows that the property (|41|) holds for k = 0.) 

We note that v G C'^+^(R), because v' = w C'^+-^(]R). Moreover, we have Wn+iv^^^ = 

we also have: 

{Wn+l v'+'Y = K+1 v'^' + {k + l) {Wn+1 v") W G C''+\^) . 

Therefore, we get Wn+iv^^^ G C^^'^iM). 

Using the relation (|38|) . we immediately obtain w G 

We have w^^+i v^'^'^ = {w'^+i v^^^) v G C'^+^(]R), and using (jlOJ, we derive: 

= {dign+l){v,Wn+l) v''^^ W + {d2gn+l){v , Wn+l) w'^+l v'''^^ + {k + 1) gn+l{v , Wn+l) W , (42) 

where dign+i (resp. d2gn+i) denotes the partial derivative of gn+i with respect to its first (resp. 
second) variable. From the definition we see that gn+i{v,Wn+i) can be decomposed as 

follows: 

gn+l{v,Wn+l) = -{n + '2)wl^i v"''^^ + Wn+l Pi{v) + Po{v) , 

where Pq, and Pi are polynomial functions. Using this decomposition, and the induction as- 
sumption H4ip . we can show that each term of the sum in the right-hand side of ()42p belongs 
to (7*^+^(1^). Consequently w'^+i v'''^'^ belongs to C''+^(M), and (jUJ) holds with k replaced by 
k + 1. Because (|41|) holds for k = 0, we get that (|4H) holds for A; = n + 1, so we have proved 
w G C"+2(M), and v G C"+3(M). 

Step 4: it remains to show that w satisfies the asymptotic expansion (|34|) at the order n + 1. 
Using and Wn+i G C^W, we obtain: 

Wn+l{0 - Wn+l = Wn+2 *(0 + o(*(0) > aS ^ ^ 0. 

Plugging this expansion in we obtain at the order n + 1, so the proof of the induction 
is complete. □ 

Once we know that the function v belongs to C"^^(M), for a G (0, an], then v = v+{v-+v+)/2 
also belongs to C"~''^(]R), and we have already seen in the previous section that v does not vanish 
because v{^) > vj^^ > for all ^. Moreover, the components {p, u, e) of the shock profile are given 
by: 

v{i) 7-1 

so one has (p, u, e) G C^''~^(M), and the proof of Theorem|2lis complete. (Recall that the strength 
of the shock tends to zero if, and only if a = |m+ — ii_|/2 tends to zero.) 

A Formal derivation of the model 

It is worth describing how the model © can be obtained from a more complete physical 
system. The derivation we propose below remains formal - a rigorous proof being certainly 
delicate and beyond the scope of this work - and we refer to |Bue041 IGouOGl ILMH991 IMih84j 
for further details. Let us introduce the specific intensity of radiation f{t, x, v), that depends on 
a time variable i > 0, a space variable x G M^, and a direction v G We make the 'grey 

assumption', which means that the frequency dependence is ignored (all photons have the same 
frequency). Photons are subject to two main interaction phenomena: 



scattering produces changes in the direction of the photons, 



• absorption/emission where photons are lost/produced through a transfer mecanism with 
the surrounding gas. 

The scattering phenomenon is described by the operator: 



Qs{f)it,X,v) 



gJV-l 



f(t,x,v')dv' - f{t,x,v) 



(with dv the normahzed Lebesgue measure on and the absorption/emission phenomenon 

is described by the operator: 



Qa{f ){t, X, V) = Ua e{t, X)^ - f{t, X, v)) , 



where we used the Stefan-Boltzmann emission law, 9 being the temperature of the gas, and 
cr the Stefan-Boltzmann constant. In these definitions, the coefficients as,a are given positive 
quantities. These phenomena are both characterized by a typical mean free path, denoted £s,(^a 
respectively. Therefore, the evolution of the specific intensity is driven by: 



-dtf + vVJ = y Qsif) + Y QaU) = Qif) 
C ta 



(43) 



where c stands for the speed of light. The equation (|43j) is coupled to the Euler system describing 
the evolution of the fluid: 



dt{pu) + • (pn n) + V^^P 
dt{pE)+V^ ■ {pEu + Pu) = - 



1 



Q{f)dv. 



vQ{f)dv, 



(44) 



The equations (|43|) . (|44|) are thus coupled by the exchanges of both momentum and energy, and 
by the Stefan-Boltzmann emission law. Observe that only the emission/absorption operator 
enters into the energy equation since the scattering operator is conservative (this would be 
different if Doppler corrections were taken into account) . Note also that the total energy: 



- / / f dv dx + [ pE dx , 

c Jrn JsN-1 JrJV 



is (formally) conserved. Writing the system H44() . and the kinetic equation (|43() in the dimen- 
sionless form, we can make four dimensionless parameters appear: 

- C, the ratio of the speed of light over the typical sound speed of the gas, 

- Cs, the Knudsen number associated to the scattering, 

- Ca, the Knudsen number associated to the absorption/emission, 

- V, which compares the typical energy of radiation and the typical energy of the gas. 
We thus obtain the rescaled equations: 



'Utf + V-VJ = ^ Qsif) + ^ Qaif) 

dtP + ^x ■ (pu) =0, 

V 

dt{pu) + Vx • (pu (g) u) + VxP = 7-0-. 
dt{pE)+V^ ■ [pEu + Pu) = -^<Ja'^^ 

'-a 



vf{v) dv , 
f{v) dv 



(45) 



System Q, Q is then obtained in two steps. First of all, we assume C ^ 1. Next, we keep V of 
order 1, and we are concerned here with a regime where scattering is the leading phenomenon: 
the mean free paths are rescaled according to: 



1 

c 



Ca ^ C 



The asymptotics can be readily understood by means of the Hilbert expansion: 

Identifying the terms arising with the same power of 1/C, we get: 

- at the leading order, f^^^ belongs to the kernel of the scattering operator, so that is does not 
depend on the microscopic variable v: f^^\t,x,v) = n{t,x), 



the relation Qs(/^^^) = v • V^/(°) then leads to: f'^'^\t,x,v) 

integrating the equation for /^^^ over the sphere yields: 

1 



dtn 



A^n = (Ja{6 -n) 



Note also that in the momentum equation, we have: 

^/ vf{v)dvc:^as 1 V f^'^\v)dv = -^VxTi. 
Finally, we obtain the limit system: 



dtp + Vx- (pu) =0, 



V 



dtipu) + Vx ■ {pu® u) + V^P = -j^ ^xfi, 
dt{pE)+Vx-{pEu + Pu) = -Vaa {t 



n) 



(46) 



dtn 



n] 



The system (|46() describes a nonequilibrium regime, where the material and the radiations have 
different temperatures {9 ^ n^^^); the equilibrium regime would correspond to assuming that 
the emission/absorption is the leading contribution. 

After this first asymptotics, we perform a second asymptotics where we set: 

Vaa = l, Nas = l/aa. 

This leads to ©. Of course, one might wonder how this second approximation modifies the 
shock profiles compared to H46|) . in particular when we get rid of the radiative pressure in the 
momentum equation. We refer to |Mih84l page 579] for some aspects of this problem. 
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